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Ab initio calculations of edge-functionalized armchair graphene nanoribbons: Structural, 

electronic, and vibrational effects 

Nils Rosenkranz/'|j Christian Till/ Christian Thomsen/ and Janina Maultzsch^ 
^ Institut filr Festkorperphysik, Technische Universitdt Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 

We present a theoretical study on narrow armchair graphene nanoribbons (AGNRs) with hydroxyl function- 
alized edges. Although this kind of passivation strongly affects the structure of the ribbon, a high degree of 
edge functionalization proves to be particularly stable. An important consequence of the geometric deviations 
is a severe reduction of the band-gap of the investigated 7-AGNR. This shift follows a linear dependence on 
the number of added hydroxyl groups per unit cell and thus offers the prospect of a tunable band-gap by edge 
functionalization. We furthermore cover the behavior of characteristic phonons for the ribbon itself as well as 
fingerprint modes of the hydroxyl groups. A large down-shift of prominent Raman active modes allows the 
experimental determination of the degree of edge functionalization. 



I. INTRODUCTION 

The unique physical properties of graphene have aroused 
steadily growing research efforts since 2004 equally in basic 
research and applications engineering [1]. Despite its excep- 
tional electronic characteristics, the gap-less band structure of 
graphene impedes a direct use in many potential nanoelec- 
tronic devices. A possible way to overcome this obstacle is 
to design graphene nanoribbons (GNRs), which open a gap 
due to the lateral confinement of the electronic wave function 
[2,3]. The fundamental properties of these narrow stripes of a 
graphene sheet strongly depend on the edge structure as well 
as the ribbon width [3-5]. GNRs are therefore classified into 
armchair-edge (AGNR) and zigzag-edge ribbons (ZGNR) in 
analogy to the achiral types of carbon nanotubes (CNTs). Var- 
ious approaches to produce GNRs have been suggested rang- 
ing from lithographic methods to the longitudinal unzipping 
of CNTs [6-8]. However GNRs prepared in these ways have 
in common one disadvantage well-known from the production 
of CNTs: It is hard to extract a sample of defined edge struc- 
ture and width, similar to the typical chirality distributions in 
CNT outputs. In 2010, Cai et ah succeeded in synthesiz- 
ing one specific type of GNR with defect-free edges and a 
well-defined width [9]. The preparation of this 7-AGNR rep- 
resents a major breakthrough in the field as the availability of 
predefined ribbons is crucial in order to make use of many 
attributes of these compounds. However, it should be noted 
that the reported synthesis is an exception as other GNRs so 
far cannot be prepared in a comparably accurate way. Typical 
lithographic fabrication includes treatment with e. g. oxygen 
plasma which makes it very likely that GNRs are passivated 
other than with hydrogen. Diverse functional groups may be 
a side product of the manufacturing process. In this case it 
is vital to know how the physical properties of the GNR are 
affected. On the other hand, functionalization might be intro- 
duced intentionally in order to tune certain properties. Since 
cutting a graphene sheet into ribbons leaves dangling bonds, 
the edges hold a great potential for various chemical mod- 
ifications. In order to describe pristine GNRs, the saturation 
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with hydrogen is commonly assumed throughout the literature 
whereas divers configurations have been shown to be stable 
[10]. Beyond that other types of edge passivation open up an 
enormous number of ways to influence the attributes of GNRs. 
Previous works on this topic focus on ZGNRs, which yield 
edge states near the Fermi-level when passivated with hydro- 
gen [11-14]. The electronic structure of ZGNRs is thus very 
sensitive to edge modifications. To the best of our knowledge, 
edge functionalization in AGNRs has only been considered 
as a side remark by Cervantes-Sodi et al who do not expect 
major impacts due to missing impurity levels in the band-gap 
[11]. However our results suggest that also other mechanisms 
like functionalization-induced strain play an important role in 
this regard. 

In this paper, we discuss the effects of edge passivation with 
hydroxyl groups on structural, electronic, and vibrational 
properties of 7-AGNRs by means of ab initio calculations. 
We demonstrate that passivation of the ribbon edges with OH- 
groups is energetically favorable over H passivation, with a 
high degree of functionalization being likely (Sec. |llll). We 
further present deviations in the geometries of functionalized 
ribbons and their hydrogen terminated counterparts. Section 
HVl deals with the effect of hydroxyl functionalization on the 
electronic band structure of 7-AGNRs. We find a strong de- 
crease of the band-gap with increasing degree of functional- 
ization, which is of great interest with regard to possible appli- 
cations. The influence of functional groups on characteristic 
vibrational modes is described in Sec. |Vl Besides the char- 
acteristic Raman active phonons of GNRs we also consider 
modes specific to the hydroxyl groups. 



II. COMPUTATIONAL DETAILS 

The 7-AGNR is an ideal model system to explore the effect 
of functionalization not only due to its availability in perfect 
quality. It also represents an appropriate compromise between 
growing computational cost for wide ribbons and vanishing 
ribbon-like characteristics for too narrow ones. We simulated 
7-AGNRs with a varying number of hydroxyl groups at the 
edges, whereas the rest of the edge atoms was passivated with 
hydrogen. As we discuss exclusively chemical modifications 
at the edges, we define the degree of functionalization (DOF) 
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FIG. 1. (a) Unit cell of a 7-AGNR with one hydroxyl passivated 
edge atom. Carbon atoms are shown in black, oxygen atoms in red 
(gray), and hydrogen atoms in light gray. The blue (gray) arrows 
indicate the lattice constant c and the ribbon width w. (b) Super- 
cell of a 7-AGNR allowing more diversified configurations of edge 
functionalization. The example shown here has a degree of edge 
functionalization of 5/8. 



as the number of hydroxyl groups divided by the number of 
available carbon edge atoms. For one unit cell the DOF can 
take integer multiples of 0.25 (see Fig. [T]). In order to de- 
scribe intermediate values we also performed calculations of 
a super-cell containing two unit cells along the ribbon axis. 
All results presented in this work were derived by means of the 
density functional theory package SIESTA in the local-density 
approximation (LDA) [15,16]. This code does not treat the 
core electrons explicitly but within norm-conserving pseu- 
dopotentials. Furthermore strictly localized basis sets were 
used to describe the valence electrons. We chose a double- 
C, basis set the confinement of which was given by a maxi- 
mum energy shift of 30 meV. These two characteristics make 
the SIESTA method particularly efficient for low-dimensional 
systems like GNRs. In order to prevent interaction between 
neighboring images of the ribbon, the lattice vectors were cho- 
sen such that the distance between two ribbons was at least 
14 A exceeding twice the maximum radii of the basis wave 
functions for all elements. The integration grid in real space 
is ruled by an internal SIESTA parameter which was chosen 
fine enough to represent plane electron waves of up to 350 Ry. 
Since continuous states only appear along the ribbon axis we 
fixed the sampling grid of the Brillouin zone to 1 x 1 x 50. 
We relaxed the coordinates of all structures until the maxi- 
mum atomic force was less than 0.01 A. The frozen phonon 
approximation was used to calculate the force constant ma- 
trix. Following the procedure in previous works we scaled 
the r point frequencies of all modes by a constant factor in 
order to compensate the underestimation of the bond lengths 
in LDA calculations [5,17,18]. The scaling factor of 0.974 
found for the E2g mode of graphene matches the discrepancy 
between LDA simulations and Raman measurements of the 
high energy band of 7-AGNRs [9]. 



III. STRUCTURE AND STABILITY 

The addition of hydroxyl groups to the edges of 7-AGNRs 
induces considerable geometric effects on the ribbon itself as 
discussed in the following. To get a clear picture we sim- 
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FIG. 2. Out-of-plane arrangement of the functional groups in a 7- 
AGNR. The bottom picture shows the view along the ribbon axis as 
indicated by the black arrow. 



ulated a variety of configurations with different degrees of 
functionalization. Apart from calculating a single unit cell 
of a 7-AGNR we also considered a super-cell containing two 
unit cells as shown for two examples in Fig. [U We find 
the planar configuration of AGNR and hydroxyl group to be 
metastable in contrast to findings on hydroxyl functionalized 
ZGNRs [11]. A bending of the functional groups out of the 
ribbon plane is energetically favorable by 0.12 - 0.14eV per 
group, depending on the degree of functionalization. For a 
single unit cell (Fig. [1^) our results show an opposite dis- 
placement of the functional group and its neighboring edge 
site regardless whether the latter one is passivated with hydro- 
gen or another OH group (see Fig. O. In case of the 1 x 1 x 2 
super-cell, the lower translational symmetry leads to an oppo- 
site bending of entire carbon hexagons at one edge if at least 
one functional group is involved. Apart from the out-of-plane 
arrangement, the edge passivation with OH groups causes a 
stretching along the ribbon axis accompanied by a squeezing 
across the width. The relative variations of the lattice con- 
stant c and the ribbon width w are shown in Fig.O Obviously 
the strain on the unit cell increases considerably with growing 
degree of functionalization. Although the total width of the 
ribbon is strongly decreased, the C-C bonds along the width 
are shortened only close to the edges. In contrast, central C-C 
bonds in this direction are even longer in the functionalized 
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FIG. 3. Relative variation of the lattice constant (black triangles) and 
the ribbon width (red squares) for different degrees of functionaliza- 
tion. Solid lines are guides to the eye. 
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FIG. 4. Binding energy per hydroxyl group depending on the de- 
gree of functionalization. Red diagonal crosses represent data de- 
rived from calculations of one unit cell. Black crosses are values 
found for a super-cell containing two unit cells along the ribbon axis. 
The solid black line is a guide to the eye. 



case than in the pristine ribbon. It turns out that the overall 
deformation of the unit cell of functionalized ribbons is due 
to a decrease of the bond angles along the width by about 
2°. Altogether, we observe a remarkably high Poisson ra- 
0.65. We performed additional calculations 
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on non-functionalized 7-AGNRs under uniaxial strain which 
yield a much lower Poisson ratio of 0.25 [18]. This distinction 
has to be taken into account when comparing functionalized 
and strained GNRs as discussed in Sec. [IVl 
Since the chemical stability of the given structures is crucial 
for any experimental realization, we calculated the binding en- 
ergy per OH group Eb for all investigated configurations as 
follows: 
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£^(GNR/) and E{GNRH-term) are the total energies of 
a ribbon passivated with Nf hydroxyl groups and a fully H- 
terminated ribbon, respectively. £"(0) represents the total en- 
ergy of an isolated oxygen atom. The resulting binding en- 
ergies from calculations of both the super-cell and unit cell 
approaches are shown in Fig. ID First of all we observe a sub- 
stantial energy gain for the edge hydroxy lation of 7-AGNRs in 
agreement with studies of the same effect in ZGNRs [12]. The 
binding energy clearly increases with growing degree of func- 
tionalization. Therefore a fully functionalized ribbon repre- 
sents the most stable configuration as found for ZGNRs [12]. 
In many cases it is possible to construct inequivalent configu- 
rations with the same number of OH groups which then yield 
differing energy values. This effect is most pronounced in 
the unit cell approach for a degree of functionalization of 0.5. 
All studied cases show a clear energetic preference of single 
edge functionalization. Particularly high stabilities of fully 
saturated edges indicate the occurrence of hydrogen bonds be- 



tween neighboring OH groups which have also been predicted 
for hydroxy lized ZGNRs by Hod et al. [12]. Figured further 
shows consistently higher energy gains for super-cell simu- 
lations. This agrees well with the above discussed enhanced 
possibilities for an out-of-plane arrangement - and thus for the 
energy minimization - at lower translational symmetry restric- 
tions. 



IV. ELECTRONIC BAND STRUCTURE 

A high potential of various edge modifications for altering 
the band structure of ZGNRs has been reported previously 
[11-13]. Although AGNRs do not have edge states near the 
Fermi level in contrast to ZGNRs our studies reveal drastic 
modifications of the band structure of 7-AGNRs under edge 
functionalization with OH groups. Figure [5] shows the elec- 
tronic bands of a pristine and a maximally functionalized 7- 
AGNR. The most striking result is the shift of the bands close 
to the Fermi level which leads to a strong decrease of the band- 
gap. We find a linear dependence of the band-gap on the de- 
gree of functionalization as shown in Fig. [6l Inequivalent 
configurations with the same degree of functionalization give 
different values in analogy to the discussion of the total energy 
in Sec. Hill Furthermore, simulations of one unit cell yield 
slightly lower band-gaps than super-cell calculations. Never- 
theless both approaches result in a perfectly matching linear 
shift of the band-gap with increasing degree of functionaliza- 
tion. Altogether the band-gap of 7-AGNRs can be reduced by 
~ 0.7eV or almost 50 % by means of edge functionalization 
with hydroxyl groups. This is highly interesting with regard 
to nano-electronic applications as it offers in principle a tun- 
able band-gap over a wide range. We did not take into account 
electron-electron and electron-hole interaction [19,20]. How- 
ever, we believe that the general behavior will be conserved, 
as the dominant contribution to the band-gap shift is from ge- 
ometrical effects as discussed in the following. 

The observed band-gap dependence very much resembles 
previous results on uniaxially strained AGNRs, suggesting an 
interpretation based on the geometric deformation of the rib- 
bon produced by the hydroxylation [18]. However, tensile 
strain along the ribbon axis of 2 %, which is the maximum 
deviation of the lattice constant under functionalization, corre- 
sponds to a decrease of only 0.2 eV in band-gap. Considering 
the extraordinarily high Poisson ratio caused by the OH passi- 
vation of the ribbon edges, we directly compare the variations 
in width and band-gap. The width reduction of 1 .4 % observed 
for the fully functionalized 7-AGNR (see Fig. [3]) would cor- 
respond to a tensile strain along the ribbon axis of 5.6 %. This 
high strain would give rise to a band-gap reduction of 0.6 eV 
which comes very close to the effect found for functionalized 
ribbons. Based on this analysis and on the fact that AGNRs of 
other widths show an opposite strain dependence of the band- 
gap, we assume that other AGNRs show a contrary band-gap 
behavior under functionalization with hydroxyl groups [18]. 
Given the successful preparation of smaller or wider AGNRs, 
this would additionally offer a material with a tunable band- 
gap in opposite directions. 
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FIG. 5. Band structures of an H-terminated 7-AGNR (solid black 
line) and a 7-AGNR with every edge atom being hydroxyl passivated 
(dashed blue line). 
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width- sensitive but independent of the edge type. Edge- 
functionalization distorts the regular displacement pattern of 
the BLM considerably in the non- symmetric cases with one 
or three OH groups per unit cell. Nevertheless the basic 
breathing-like pattern is preserved in all investigated exam- 
ples. Figure |7ja) shows the BLM frequency of 7-AGNRs de- 
pending on the degree of functionalization. Obviously the vi- 
bration is substantially damped by the OH groups, resulting in 
a red-shift of 95 cm~^ in total. As the BLM is Raman active, 
this huge down- shift allows to determine the degree of func- 
tionalization experimentally. However in practice this goal 
may be hindered by the plenitude of Raman active modes in 
less symmetric compounds. The great influence of the partic- 
ular arrangement of the hydroxyl groups on the stability and 
the band-gap is also reflected in the energy of the BLM (see 
Fig- Eta) at a degree of functionalization of 0.5). In previ- 
ous work we showed that uniaxial strain also provokes a lin- 
ear red-shift of the BLM. However, the dimension of the ob- 
served shift under edge passivation does not allow an exclu- 
sive ascription to the functionalization-induced strain even if 
we consider the unusually large Poisson ratio as described in 
Sec. HVl Most likely multiple factors contribute to the strong 
frequency shift. Thus, the addition of functional groups to 
the ribbon edges may produce a similar effect as a broadening 
of the ribbon which is known to reduce the BLM frequency. 
This assumption is supported by the displacement pattern of 
the BLM shown in Fig. |7tb). All atoms of the functional 
groups largely follow the breathing-like motion of the ribbon 
itself and vibrate in phase with the carbon atoms perpendic- 
ular to the ribbon axis and in the ribbon plane. Within the 
picture of an expanded ribbon structure, the oxygen atoms 
stand for an additional dimer at each edge of the ribbon. Thus 



FIG. 6. Band-gap of hydroxyl-passivated 7-AGNRs depending on 
the degree of functionalization. Red (gray) vertical crosses and black 
crosses represent band-gap energies found for calculations of a unit 
cell and for a 1 x 1 x 2 super-cell, respectively. 



V. VIBRATIONAL PROPERTIES 

In the following we investigate the behavior of the most 
important vibrational modes of 7-AGNRs under different de- 
grees of hydroxylation. We present data on the breathing-like 
mode (BLM) and the fundamentals of the transversal and lon- 
gitudinal optical modes at ~ 1600 cm~^ which dominate the 
Raman spectrum [9,17,21]. In addition, characteristic modes 
of the hydroxyl groups and the influence of the functionaliza- 
tion on the stretching modes of the remaining C-H bonds are 
discussed. 
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A. Breathing-like mode 

Similar to the radial breathing mode in carbon nanotubes 
the BLM is likely to play a key role in spectroscopic charac- 
terization of GNRs. This fully symmetric vibration is highly 



FIG. 7. (a) Down- shift of the breathing-like mode of 7-AGNRs un- 
der functionalization with hydroxyl groups, (b) Eigenvector of the 
BLM of a fully hydroxylized 7-AGNR. Carbon, oxygen, and hydro- 
gen atoms are shown in black, red (gray), and light gray, respectively. 
Blue arrows indicate the atomic displacements. 



the fully functionalized 7-AGNR can be thought of as a 9- 
AGNR. The breathing-like mode of the hydrogen-terminated 
9-AGNR is found at 300 cm~^. Given the rough underlying 
approximation this value is fairly close to the BLM frequency 
of 285 cm"^ found in the fully hydroxylated 7-AGNR. 



B. High-energy band 



stretching modes are found at 800-1 150cm~^ [22,23]. The 
coupling of the TO to the low-energetic C-0 stretching mode 
could potentially give rise to the observed softening of the TO 
frequency. In contrast, the C-O-H bending which couples to 
the LO is expected to be closer in frequency to the LO-TO re- 
gion. Therefore the LO should be affected less than the TO in 
agreement with our findings. 



The degeneracy of the G band in graphene is lifted in 
GNRs, which show the longitudinal optical phonon (LO) at 
lower energies than the transverse mode (TO) [5]. Both 
phonons are clearly measurable and distinguishable in Ra- 
man experiments though under different polarization condi- 
tions [17]. As presented in Fig. Oa) we find a fundamen- 
tally differing behavior of both modes under functionalization. 
Whereas the TO mode shifts to lower energies with increasing 
degree of functionalization similar to the BLM, the LO does 
not show such a trend. Aside from the strongly varying fre- 
quencies of the three inequivalent configurations with two OH 
groups per unit cell, the energy of the LO seems to be indepen- 
dent of the edge passivation. The large shift of the TO, in turn, 
is well- suited for a determination of the degree of functional- 
ization. Like in case of the BLM a strain-based interpretation 
of the LO and TO behavior is problematic. As shown previ- 
ously, tensile strain along the ribbon axis provokes a down- 
shift of both the LO and TO modes [18]. In contrast to our re- 
sults on functionalized ribbons, uniaxial strain affects the LO 
stronger than the TO. Furthermore, even the strain-induced 
shift of the TO exceeds the effect of hydroxylation by a factor 
of 2 [18]. Once again a look at the TO and LO eigenvectors 
shown in Fig. [8tb) is helpful. A noticeable difference be- 
tween the two displacement patterns concerns the behavior of 
the functional groups. Whereas the oscillation of the carbon 
atoms in the TO provokes a stretching of the C-0 bonds, a 
strong bending of the OH group dominates the LO pattern. In 
both cases the motion of the outermost carbon atoms induces 
these characteristic vibrations of hydroxyl groups which are 
usually found at lower frequencies. Depending on the com- 
pound which is hydroxy lized, in-plane deformations of the O- 
H groups are usually localized at 1260-1 440 cm~^ and C-0 
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FIG. 8. (a) Frequencies of the fundamentals of the LO (left) and TO 
phonons (right) of hydroxy lized 7-AGNRs. (b) Displacement pattern 
of these modes. The same color code as in Fig. [Tib) is used. 



C. Functionalization-specific modes 

Finally we discuss the high-energy range of the vibrational 
spectrum comprising the stretching modes of C-H and 0-H 
bonds. A hydrogen- terminated 7-AGNR yields two doubly 
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FIG. 9. Stretching frequencies of (a) the C-H bonds and (c) the 
0-H bonds in 7-AGNRs depending on the degree of hydroxylation. 
Vibrational frequencies found for single-edge functionalization and 
the two different cases with half-functionalized edges are represented 
with the symbols stated in b). (b) Three inequivalent configurations 
with two OH groups per unit cell. Different elements are again dis- 
played in the color code of Fig. [3b). 



degenerate pairs of C-H stretching vibrations. This number 
is reduced in functionaHzed ribbons to the number of remain- 
ing C-H bonds. Consequently we find as many 0-H stretch- 
ing modes as involved hydroxyl groups per unit cell. As dis- 
played in Fig. [3a) C-H stretching generally occurs in two 
narrow frequency ranges around 3000 and 3020 cm~^ , respec- 
tively. Figure Oc) shows three inequivalent structures with a 
degree of functionalization of 0.5. The two cases with one 
hydroxyl group per edge maintain C2 rotational symmetry or, 
alternatively, a mirror plane. Consequently each of these con- 
figurations yields a doubly degenerate C-H stretching mode 
which is found 2997 cm~^. As expected, the non-symmetric 
single-edge hydroxy lated 7-AGNR shows two separated C-H 
stretches. Measuring the C-H stretching frequencies would 
allow some insight in how many OH-groups are bound to the 
ribbon edges. From the absence of the modes at 3020 cm~^ a 
high degree of functionalization may be deduced. 
We now turn our attention to the 0-H stretching vibrations 
presented in Fig. Ob). Apart from the fully functionaHzed 7- 
AGNR all cases show such modes around 3200 cm~^. It turns 
our that these are vibrations of half functionaHzed edges, /. e. 
of OH groups which do not have direct hydroxyl neighbors. 
As soon as a completely hydroxylized edge is involved we 
observe a severe down-shift, which is a clear indication of hy- 
drogen bonds. Opposite-phase stretching vibrations of neigh- 
boring 0-H groups are found at 2605 to 2630 cm~^, whereas 
the corresponding in-phase modes are further lowered to 2120 
to 2145 cm~^. Consequently the fully hydroxylized 7-AGNR 
does not show any OH- vibrations around 3200 cm~^. The ex- 
perimental detection of the 0-H stretching modes thus allows 
to further narrow down the degree of functionalization of a 
given sample. Together with the results presented in Sees. 



IV Al and IV Bl this suggests that an extensive characterization 
of functionaHzed 7-AGNRs is feasible via vibrational spec- 
troscopy. 



VI. CONCLUSION 

We presented a thorough investigation of the structural, 
electronic, and vibrational properties of 7-AGNRs with hy- 
droxyl functionaHzed edges. The passivation of every carbon 
edge atom with an OH group gives the most stable config- 
uration. At the same time this maximum functionalization 
has the biggest effect on the electronic band structure and on 
the main vibrational modes. The most striking result with re- 
gard to possible applications is the large shift of the band-gap 
which can be attributed to strain induced by the hydroxyla- 
tion. A tunable band-gap in the given range might be very 
useful in future nanometer- sized electronic devices. Finally 
we showed that the degree of edge functionalization may be 
determined by Raman spectroscopy. These findings are par- 
ticularly exciting in the light of the successful preparation of 
pure 7-AGNRs. An experimental verification of our results is 
thus within the realms of possibility. Moreover it will be en- 
lightening to study the impact of hydroxyl groups on AGNRs 
of different widths as well as the effect of other types of edge 
functionalization . 
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